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ABSTRACT 

The dynamics of accretion onto a schwarzschild black hole is studied us- 
ing Paczynski-Wiita pseudo newtonian potential. Steady state solution of the 
flow equations is obtained using thin disc approximation, including the effect 
of bremsstrahlung cooling in the energy equation. The topology of transonic 
solutions are got and the conditions for shock formation (Rankine-Hugoniot con- 
ditions) are checked. Shock and no-shock regions in the parameter space of 
accretion rate and specific angular momentum are identified. We motivate the 
time-dependent numerical simulation of the flow, as a candidate for explaining 
the quasi periodic oscillations observed in black hole candidates. 

Subject headings: accretion, accretion disks — black hole physics — hydrody- 
namics — methods: numerical — shock waves 

1. Introduction 

The parameter space of shock formation in adiabatic transonic accretion discs was first 
discussed in Chakrabarti(1989). Though time-dependent numerical simulation of transonic 
accretion disc in the presence of cooling was dealt in Molteni, Sponholz & Chakrabarti 
(1996, hereafter MSC96), the detailed study of parameter space was reserved for the future. 
In this paper we describe the procedure for obtaining the flow topology in the presence of 
bremsstrahlung cooling and obtain the dependence of topology and shock formation on the 
parameters, accretion rate and specific angular momentum. 
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2. Flow equations 

The equations which govern the flow are the basic conservation equations of mass, 
momentum and energy. The pseudo newtonian potential of Paczynski-Wiita (1980) is used, 
which mimics the gravitational field around a schwarzschild black hole to a sufficient accuracy. 
The effect of bremsstrahlung cooling (Lang 1980) in electron-proton plasma is included in 
the energy equation. v(v r , v z ),p and p are the velocity, density and pressure respectively. 

Continuity equation: 

f + V.( P v)=0 (1) 

Euler equation: 

! + (v . V )v^-V g (2 ) 

where 

g = i S 

the Paczynski-Wiita potential, G is the gravitational constant, M is the 
mass of the black hole, r g = is the schwarzschild radius and c is the velocity of light. 

Energy equation: 

V.( P 6v) + A-r = -^(p e ) (3) 
where e = \v 2 + U + - + g + is the specific energy, U = -, is the thermal energy, 7 is 



the adiabatic index, A = 1.43 x 10" 2 ' '^T 1 ^ 2 'gj is the expression for bremsstrahlung cooling, 
T is the temperature, m p is the mass of proton, gj is the gaunt factor and T is the heating 
term. 

Thin disc approximation: 

We use cylindrical polar coordinates (r, 0, z) for the axisymmetric flow. Because of axisym- 
metry and thinness of the disc the following assumptions are made. 

p(<j>°, z°) i.e. p is not a function of <fi and z 



v z (r°,<f) ,z ) 
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p(r) i.e. p is a function of r only 

g{r) 
e(0°) 

By the help of these assumptions it is possible to write the steady state equations in the 
following form. 

V-rf/r + jgfi^TV* dv GM . 1 2 



a 2 /v — v dr ' (r — r g ) 2 

1 dp ldv 1 



(r-r g ) 2 
(5) 



p dr v dr r 

dv 2a da a 2 dp , . „ . . 

+ — - r + — -f + ^' = 6 

ctr 7 dr pj dr 

where ' denotes the derivative with respect to r, j = 1.43 x 10~ 27 in cgs units and the 
polytropic relation p = Kp 1 is used to obtain the sound speed a = \J~§j> — These 
equations are solved using fourth-order Runge-Kutta method (Press et al. 1992). 



3. Sonic point analysis 

The flow can pass through a point where denominator in the expression for ^ becomes 
zero. If it happens that the numerator also becomes zero, then ^ can be finite. We call such 
a point a critical point. In this problem the critical point is same as sonic point. At critical 
point, using 1' Hospital's rule, we get (for 7 = |), 



dv -B ± y/(£? 2 - AAC) 

dr 2A [ ' 



A=l (8) 



B = - 1 - V --k P a 2a - 2 (l + -) + -^ (9) 
3r V 3 J 3 v K ' 
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c = \^ + r - |c— p(i - «) - t«^ 2q -v 

3 r 3 r 9 



(10) 




) 



a 



(11) 



where a = 0.5, [i = 0.5 and k is the boltzmann constant. Now starting from the critical 
point we obtain the transonic solution topologies. 



Fluid dynamical problems are inherently sensitive to the boundary conditions. We 
formulate the problem as follows. The variables which are functions of space and time are 
v,p,p or v, p, a. The parameters are accretion rate, specific angular momentum(A), 7, 
cooling process and viscosity. We consider inviscous flow, set 7 = 5/3 and the mass of the 
black hole M is chosen as 1O 8 M . So the free parameters are accretion rate and A. The 
location of the critical points, r cl ,r c2 , is the freedom we have, subjected to the constraint 
that shock conditions should be satisfied. The natural constraint of stability is likely to 
decide the uniqueness of the solution. The algorithm for obtaining the numerical solution is, 

1. Choose accretion rate and A 

2. Choose r cl and r c2 range and find corresponding a c \ and a c2 

3. Obtain solution topology by using fourth-order Runge-Kutta method 

4. Check the shock conditions 

The typical components of the topology are shown in fig.l, where mach number is plotted 
as a function of the radial distance. Fig. 2 and 3 shows the solution topologies for chosen 
parameter values. Rankine-Hugoniot shock conditions (Landau & Lifshitz 1984) basically 
ensure that mass, momentum and energy are conserved in spite of a discontinuity. For 
infinitesimally thin and non-dissipative shock, the shock conditions take the form, 



4. 



Numerical scheme 



Pi + PlVl 2 =P2+ P2Vl 2 , 




7-1 2 




(12) 



where the suffixes 1 and 2 refer to pre-shock and post-shock quantities. 
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5. Parameter space 

We use the following procedure to obtain the parameter space. The accretion rate in 
eddington units, is varied in the range (.0001, 500) and A in the units of is varied 

in the range (1.6,2.5). It is suggested that (Chakrabarti 1996) flows with A in the range 
between that of marginally stable and marginally bound orbit would form steady accretion 
discs, when self-gravity of the disc is neglected. For a chosen accretion rate and A we scan 
the r-axis from 1.5r 9 to 1000r 9 to find if it can be a critical point. The critical point can 
be of X or Alp or plA or 1A or x type as shown in fig.l. X looks like the English alphabet 
X, Alp like Greek alphabet alpha which opens towards infinity, pi A is reflected Alp which 
opens towards inner boundary, 1A is plA which doesn't reach the inner boundary when A is 
high and x is X which doesn't reach the inner boundary when accretion rate is high. For 
certain parameter values there is an intermediate range of r which cannot be a critical point 
(denoted by - in topology column of Table 1). Table 1 shows the parameter space. We 
obtain supersonic branches for outer critical points and subsonic branches for inner critical 
points. When shock conditions are satisfied the flow makes a transition from supersonic 
to subsonic branch and reaches the inner boundary (chosen as 1.5r 9 ) supersonically. When 
shock conditions are not satisfied accretion is still possible if X type critical point exists. 



6. Discussion and conclusions 

Of all the different possible branches of accretion for chosen parameter value, the real 
flow is likely to choose the branch which is most stable, as perturbations are always present in 
a real situation. If the assumption of steady flow is relaxed the flow might still choose a steady 
solution branch if the time scale of change is large. Time-dependent numerical simulation of 
the flow (MSC96) shows that flow solution oscillates about the steady state solution when 
certain resonance condition is met. The 'perturbations' which occur in a numerical code and 
physical perturbations should be related, to increase the faith in numerical results. Such 
numerical studies will be pursued in the future. The oscillation of shock location (MSC96) 
would mean the size of post-shock region, which is the source for hard photons (Chakrabarti 
& Titarchuk 1995), is also oscillating. These would result in quasi periodic oscillations as 
reported in Chakrabarti & Manickam (2000). 

We thank Indian Space Research Organization for funding the project, Quasi-Periodic 
Oscillations in Black Hole Candidates, of which this work forms a part. 
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Table 1. 



accretion rate 


A 


topology 


shock 


.0001 


1.6 


X 




.0001 


1.65 


T 7" 1 A A 1 T 7" 

X plA - Alp X 


noshock 


.0001 


1.7 


X plA - Alp X 


shk 


.0001 


1.9 


"\7" 1 A A 1 "\7" 

X plA - Alp X 


shk2 


.0001 


2.0 


X - Alp 


noshock 


.001 


1.6 


X 




.001 


1.65 


X plA - Alp X 


shk 


.001 


1.9 


AT" 1 A A 1 "\7" 

X plA - Alp X 


shk2 


.001 


2.0 


X - Alp 


noshock 


.01 


1.6 


X 




.01 


1.65 


X plA - Alp X 


shk 


.01 


1.9 


X plA - Alp X 


shk2 


.01 


2.0 


X 1A - Alp 


noshock 


.1 


1.6 


X 




.1 


1.65 


X plA - Alp X 


shk 


.1 


1.9 


X plA - Alp X 


shk 


.1 


2.0 


X - Alp 


noshock 


1.0 


1.6 


X 




1.0 


1.65 


X plA - Alp X 


shk 


1.0 


1.7 


X plA - Alp X 




1.0 


1.8 


X plA - Alp X 




1.0 


1.9 


X - Alp X 


shk2 


1.0 


2.0 


X - Alp 


noshock 


1.0 


2.1 


X - Alp 




1.0 


2.3 


X 1A - Alp 




1.0 


2.4 


1A - Alp 




1.0 


2.5 


1A - Alp 




2.0 


1.6 


X 




2.0 


1.65 


X plA - Alp X 


shk 


2.0 


1.8 


X plA - Alp X 




2.0 


1.9 


X 1A - Alp X 


shkl 



- 11 - 







Table 1- 


-Continued 


a pprp t, ion rat.p 


\ 


tonnlnp'v 

Lu i juiu &,y 


shock 


2.0 


2.0 


X 1A - Alp 


noshock 


2.0 


2.3 


X 1A - Alp 




2.0 


2.4 


x 1A - Alp 




5.0 


1.6 


X 




5.0 


1.65 


X plA - Alp X x 


shk 


5.0 


1.7 


X plA - Alp X 




5.0 


1.8 


X plA - Alp X 




5.0 


1.9 


X 1A - Alp X x 


shk 


5.0 


2.0 


X 1A - Alp X x 


noshock 


5.0 


2.3 


X x 1A - Alp x 




20. 


1.65 


X plA - Alp X x 




20. 


2.4 


X x 1A - Alp x 




30. 


1.65 


X plA - Alp X x 




50. 


1.65 


X x 




50. 


1.7 


X plA - Alp X x 


shk 


50. 


1.8 


X plA - Alp X x 


shkl 


50. 


1.9 


X plA 1A - Alp X x 


shk 


50. 


2.0 


X x 1A - Alp X x 


noshock 


50. 


2.4 


X x 1A - Alp x 




500. 


1.6 


X x 




500. 


1.8 


X x 




500. 


2.4 


X 





